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We present a collection of new, open-source computational tools for numerically modeling recent 
large-scale observational data sets using modern cosmology theory. Specifically, these tools will allow 
both students and researchers to constrain the parameter values in competitive cosmological models, 
thereby discovering both the accelerated expansion of the universe and its composition (e.g., dark 
matter and dark energy) . These programs have several features to help the non-cosmologist build an 
understanding of cosmological models and their relation to observational data: a built-in collection of 
several real obervational data sets; sliders to vary the values of the parameters that define different 
cosmological models; real-time plotting of simulated data; and calculations of the goodness 
of fit for each choice of parameters (theory) and observational data (experiment). The current 
list of built-in observations includes several recent supernovae Type la surveys, baryon acoustic 
oscillations, the cosmic microwave background radiation, gamma-ray bursts, and measurements of 
the Hubble parameter. In this article, we discuss specific results for testing cosmological models 
using these observational data. These programs can be found at http://www.compadre.org/osp/ 
Items/detail . cf m?ID=12406| 

PACS numbers: Pacs here 



I. INTRODUCTION 



In recent decades, the field of cosmology — both ob- 
servational data and theoretical models — has provided 
two very significant insights regarding the nature of our 
universe. One discovery, which earned the 2011 Nobel 
Prize in Physici^lf^, is that our universe is not just ex- 
panding; it is expanding at a rate that is increasing with 
time. (I.e., the matter in the universe is experiencing 
a repulsion that overcomes the attractive force of grav- 
ity.) The other major discovery is that — in order to 
successfully model the (numerous) recent astronomical 
observations — our universe must be composed of mostly 
dark matter and dark energy, with only 4% ordinary mat- 
ter (e.g., atoms Clearly, these results are of intrinsic 
interest and should be understood by people outside the 
field of cosmology. In fact, several recent articles in non- 
specialist journals have discussed these latest findingJ^M^. 
Unfortunately, resources have not existed that allow the 
broader physics community (non-cosmologists) to appre- 
ciate how cosmological observations inform cosmological 
models, ultimately leading to these new insights. In the 
present work, we seek to address this need. 

In order to help non-cosmologists to understand how 
recent observations lead to the discoveries described 
above, we provide a medley of open-source, user-friendly 
cosmological modeling programs, which we will refer to 
as "CosmoEJS"^. These programs allow the user to im- 
mediately become an amatuer cosmologist by fitting the- 
oretical models to the actual experimental data and vi- 
sually observing how well each model agrees with the 
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various observational data. Key features of CosmoEJS 
include a built-in collection of several real obervational 
data sets; sliders to vary the values of the parameters that 
define different cosmological models; real-time plotting of 
simulated data; and calculations of the goodness of fit 
for each choice of parameters (theory) and observational 
data (experiment). Taking advantage of these features, 
CosmoEJS has already been used with a variety of non- 
cosmologist audiences to bridge the gap between modern 
cosmology and mainstream physics. 

In Sec. [n] we describe some specific examples of mod- 
eling cosmology with experimental observations using 
CosmoEJS. We provide a brief introduction to cosmol- 



ogy in Sec. Ill and a more detailed introduction is pro- 
vided as an online appendix in Ref. j3, The necessary 
mathematical quantities and equations are introduced in 



Sec. Ill A and the relevant experimental observations are 
introduced in Sec. |IIIB| Sec. |IV| contains a summary of 
our results, and some details of the computer implemen- 
tation are described in Appendices \K\ and [B) 



II. MODELING COSMOLOGY 

CosmoEJS allows non-cosmologists to simulate the ex- 
pansion of the universe using various models and to com- 
pare the simulated results to experimental observations. 
Recent articles in non-specialized physics journals have 
discussed some of the complexities associated with mod- 
eling cosmology*^ ^, but these articles only considered 
specific scenarios from the small subset of cosmological 
models that permit exact analytical solutions. Moreover, 
it would be very difficult for a non-cosmologist to recreate 
or extend those results. Some web-based^ and mobile 
device applicationjil' do provide "cosmology calculators" 
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for calculating times and distances of simple models, but 
do not compare to data, or allow for the diversity of mod- 
els contained CosmoEJS. Very powerful numerical simu- 
lations for testing cosmological models and con strai ning 
the values of model parameters do already exist^^Hl^ but 
these tools have a steep learning curve, making them im- 
practical for use by non-specialists. CosmoEJS addresses 
all of these needs. It provides direct comparisons be- 
tween theory and experiment in the form of both plots 
and numbers, and it accurately carries out the complex 
mathematics without requiring technical expertise from 
the user. This allows the user to focus on developing a 
high-level understanding of cosmology, without technical 
(mathematical and computational) distractions. 

The process of using CosmoEJS is very straightfor- 
ward, and we encourage the user to download the pro- 
gram from Rcf. 8 in order to recreate and modify the 
plots that are discussed later in this section. Using Cos- 
moEJS consists of five steps: (1) loading observational 
data, (2) selecting values for the model parameters, (3) 
calculating and plotting theorectical observables, (4) as- 
sessing the goodness of fit both visually and numerically, 
and (5) plotting the expansion of the universe for the 
user-defined model. For step (1), the user can select up 
to 18 different experimental datasets, and a drop-down 
menu is provided to simplify the process of loading data. 
(See Sec. IIIB| for a description of these experimental 
data.) For step (2), users can use sliders to adjust pa- 
rameter values, subsequently changing from one model 
to another. For step (3), the calculations are carried out 
using Romberg's method of approximating integrals, and 
the user can easily vary the number of partitions to test 
for convergence of the numerical integration. For step 
(4), each time that the user changes the model's param- 
eters, several plots are generated; and for each plot, is 
automatically calculated to provide a quantitative mea- 
sure of the goodness of fit. Finally, for step (5), the size 
of the universe can be plotted versus time to see what 
type of universe results from each set of parameters; i.e., 
is the expansion of the universe constant, or accelerating, 
or decelerating? 



A. Fitting the Model 

In this subsection, we demonstrate the modeling ca- 
pabilities of CosmoEJS by comparing experimental data 
(astronomical observations) with simulated data for a 
few specific examples that do not permit analytical solu- 
tions. The experimental data consist of measurements of 
Type la supernovae (SNe), the Hubble parameter, H{z), 
gamma ray bursts (GRB), baryon acoustic oscillations 
(BAG), and the cosmic microwave background (CMB). 
(Each of these observations are described in Sec. IIIB ) 
The simulated data in this section consist of three phys- 
ically different models, each described by a different set 
of parameter values. Specifically, in each model, the uni- 
verse is chosen to have the same expansion rate today, but 



different fractions of matter, flm, and dark energy, ^a: 
^a} = {0.01, 0.99}, {0.27, 0.73}, and {1.0, 0.0}. 
Throughout this section data are plotted versus red- 
shift. It is important to note that redshift, z, can be 
used as a measure of time.'^Sl Light from nearby objects 
experiences very little redshift (z « 0), and this light 
was also emitted very recently (in terms of cosmological 
timescales). Light from far-away objects experiences a 
larger redshift, and this light was emitted longer ago.l^ 
We take advantage of this redshift /time relationship in 
multiple ways. Given a theoretical model [Eq. ([!]) in 



Sec. Ill A , CosmoEJS uses redshift values to calculate 
and display both the age of the universe today and the 
"look-back" time, which refers to how long ago the light 
was emitted that is observed to have a certain redshift. 
Also, the experimental data are measured using redshift 
(x-axis) , which provides a means of "dating" these obser- 
vational data once a particular model has been selected. 

In Fig. [T] distance modulus, /z, is plotted versus z for 
two different types of observations, SNe and GRB; and 
these observations are compared to the three different 
models, (/x is a normalized measure of the distance to an 
observed object.) We note that all three models provide 
a reasonably good fit to some of the data in Fig. [T] but 
there are also differences between the three curves that 
are clearly visible, and if the plot were scaled in (zoomed 
in) for low redshift, we would see more differences. For a 
model with more matter (more gravity), the matter den- 
sity would use gravity to try to pull the universe together, 
subsequently slowing the expansion rate, so the objects 
in the universe would be closer together (bottom curve) . 
Conversely, a universe with too much dark energy would 
expand the universe too quickly, and objects would be at 
greater distances than what is observed (top curve). 



SNe (Blacit) : X' 
GRB (Red): X' 
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FIG. 1: (color online) CosmoEJS output showing three dif- 
ferent models (curves) and two different experimental data 
sets, SNe (black) and GRB (red). The three models dif- 
fer only in their fractional matter and dark energy densities, 
{STi™, Qa} = (0.01, 0.99}, {0.27, 0.73}, {1.0, 0.0} correspond 
to blue, green, and red, respectively. (Note: Labels are drag- 
gable in CosmoEJS.) 

In Fig. [2j the expansion rates, H(z), of different galax- 
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ies are plotted as a function of their redshift, z. These 
data are compared with the same three models that are 
used for Fig. [TJ and again, all three the models have 
been defined to have the same expansion rate today, i.e., 
i/(z) = for z = 0. From these data it is clear that 
the middle curve (r2„j — 0.27, JIa = 0.73) gives a bet- 
ter fit than the other two models. This plot correctly 
displays that a model universe comprised of mostly dark 
energy (bottom curve) would have a slightly increasing 
expansion rate as a function of redshift for the range of 
redshift seen in Fig. [2] whereas a universe with mostly 
matter would have an expansion rate that drastically de- 
creases with increasing time (decreasing z). 
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FIG. 3: (color online) CosmoEJS output showing three dif- 
ferent models (curves) compared to data sets of the BAO 
ratio, the size of the sound horizon, r^, to its effective 
distance, Z)„ in the galaxies. The three models differ 
only in their fractional matter and dark energy densities, 
{fi™, Oa} = {0.01, 0.99}, {0.27, 0.73}, {1.0, 0.0} correspond 
to blue, green, and red, respectively. (Note: Labels are drag- 
gable in CosmoEJS.) 



verse, fife, and different types of dark energy {wo,Wa} 
models. 



FIG. 2: (color online) CosmoEJS output showing three dif- 
ferent models (curves) compared to a data set of the expansion 
rates H(z) of galaxies at different redshift. The three models 
differ only in their fractional matter and dark energy densi- 
ties, {f2™, Oa} = {0.01, 0.99}, {0.27, 0.73}, {1.0, 0.0} corre- 
spond to blue, green, and red, respectively. (Note: Labels are 
draggable in CosmoEJS.) 



In Figure |3] the three models show a clear difference 
when compared to the BAO ratio. Physically, the BAO 
ratio reflects the size of the sound horizon, r,, for the 



early universe baryon decoupling (see Sec. IIIB) to its 
effective distance, in the galaxies today. In Figs, [l] 
and[2j the two extreme models (upper and lower curves) 
at least fit some of the data, but they do not come close to 
the extremely precise (small error bars) BAO ratio data. 
With the addition of this third observation, it is obvious 
that only the middle curve fits well to all of the com- 
plimentary data sets. (By "complimentary", we mean 
that the cosmological parameters must be consistent with 
different observations that constrain different theoretical 
observables.) 

The user can easily vary additional parameters, sub- 
sequently adding to the complexity of the models, by 
simply adjusting the sliders for the different parameter 
values. For example, by varying the value of the current 
expansion rate, i/o, while looking at the same composi- 
tions, the changing of the expansion rate will uniformly 
scale all of the theoretical data points. In addition, the 
user can choose different spatial curvatures for the uni- 



B. Cosmological Interpretation 

For the cosmological intepretation of these fits, we pro- 
ceed to plotting the evolution of the expansion of the 
univese versus to time, t (in Gyrs), rather than redshift, 
z. Redshift is a model-independent measurement, which 
makes it an ideal quantity for the fitting that was done 
in the previous subsection. However, redshift has a non- 
linear, model-dependent relationship with time, which 
makes it very difficuff to physically interpret data that 
are plotted versus z;-— ^ For this reason, we now plot the 
expansion of the universe versus time for each of the 
three models that were studied in Sec. |II A| Specifically, 
we plot the dimensionless ratio a{t) / a{ttoday) , where a{t) 
represents the expansion function or radius (size) of the 
universe. It is also referred to as the "scale factor", and 
mathematically, it is defined as a = 1/(1-1- z)W^ 

In Figure|4j we provide plots of all three models studied 
in Sec. [Tni These plots can be interpreted to give phys- 
ical insight to the expansion at a particular time. We 
know from the fitting in Sec. |II A| that the model with all 
matter density did not match the data, and if we look at 
its expansion, we can see the slope decreases monotoni- 
cally with increasing time, so the rate of the universe's 
expansion is decreasing, which would be caused by the 
gravitational pull of the matter dominating the universe. 
This slowed expansion yields an age of only ~ 9 Gyrs. 
For the other extreme case studied in Sec. |II A[ we see a 
very early inflection point where the expansion rate went 
from decreasing to accelerating, yielding a much older 
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FIG. 4: (color online) CosmoEJS output showing three 
different models (curves) of the expansion of the universe, 
a{€) I aittoday) (scale factor), versus time. The three models 
differ only in their fractional matter and dark energy densi- 
ties, {Q„n, ^h) = {0.01, 0.99}, {0.27, 0.73}, {1.0, 0.0} corre- 
spond to blue, green, and red, respectively. (Note: Labels are 
draggable in CosmoEJS.) 



III. COSMOLOGICAL BACKGROUND 

The cosmic acceleration of the universe can be ex- 
plained by a cosmological constant, or some other form 
of repulsive dark energy, i.e. a negative pressure and a 
negative equation of state, or by an extension or modifi- 
cation to gravity at cosmological scales of distance^^^. In 
the context of general relativity (GR), to account for this 
dark energy effect, the addition of a A term to Einstein's 
Field Equations (EFE) can be used to derive equations of 
motion with a cosmological constant of the desired value 
consistent with the dynamics of Friedmann-Lemaitre- 
Robertson- Walker (FLRW) univers^^. We provide a 
means of testing this commonly accepted model of the 
Universe and others with observations of SNe Type la, 
gamma-ray bursts (GRB), baryon acoustic oscillations 
(BAO), the distance to the last scattering surface of 
the cosmic microwave background (CMB) radiation, and 
measurements of the Hubble expansion rate parameter, 
H{z), thereby deriving the parameters for the standard 
model in cosmology. 



universe of 28 Gyrs. This acceleration is caused by the 
domination of the dark energy as the universe expands. 
Finally, for the model which fit all of the observations in 
Sec. |II A[ we see an inflection point at about ~ 9 Gyrs, 
which is caused by the domination of the dark energy 
as the universe expands and the gravitational pull from 
matter is weakened. This model calculates a current age 
of the universe to be ^ 14 Gyrs. 

According to t he c ombinations of observations in the 
current literatur^^IlHII^ l^j^g best-fit values for the pa- 
rameters in CosmoEJS are {Hq ^ 70.0 km/(s Mpc), 
fif, ~ 0.04, fie - 0.23, - 0.27, ~ 0.73, ~ 0.0, 
Wo ^ —1.0, Wa ^ 0.0}, representing the so-called A Cold 
Dark Matter (ACDM) model. With CosmoEJS, it is pos- 
sible to find these values by systematically trying dif- 
ferent sets of parameters with combinations of the data 
sets. Then, a physical interpretation of the model's fit 
throughout its evolution can be made to compare with 
the cosmological observation^^. 



A. Mathematics 

In this section, we define the mathematics behind 
the theoretical models involved in CosmoEJS. Specif- 
ically, the programs assume a big bang physical uni- 
verse, a mathematical model according to general rel- 
ativity (GR), and a uniformly distributed spacetime in 
all directiona^Sj. From these assumptions, the programs 
numerically integrate an equation of motion for the dy- 
namical evolution of the expansion rate of the universe^S. 
This equation of motion can be expressed in terms of the 
Hubble expansion rate, H{z), as a function of redshift, 
z. [The theoretical details of the integration of H{z) for 
a particular observation are described in Ref. 'W, and the 
numerical implementation is shown in Fig. |6]of the Ap- 
pendix |Aj] Explicitly, we use 



Hiz) = i/o^a„(l + + Ha [(1 + z)3(i+»o+».) cxp (^^)] + W + zY- 



(1) 



This equation contains all of the parameters that can be 
varied in CosmoEJS. (All parameters are dimensionless 
except Hq). Briefly, these parameters and their currently 
accepted values are: 

• ~ 70.0 km/(s Mpc) : the Hubble Constant; 

• ilm ~ 0.27 : the fractional matter density (subject 



to the constraint: fi^ 



• ilft « 0.04 : the fractional baryon density; 

• ilc ~ 0.23 : the fractional cold dark matter density; 

• ~ 0.73 : the fractional dark energy density; 

• rjfc « 0.0 : the fractional curvature density; 
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• fio ~ 1-0 : the sum total energy density (subject to 
the constraint: JIq = + + ^k)] 

• wq ~ —1 : the equation of state of dark energy ; 

• Wa ~ 0-0 ■ the derivative of wq; 

As the values of these parameters change, Eq. ([T]) de- 
scribes different types of evolutions for the universe. The 
details of these parameters are further explained in the 
next paragraph. 

The Hubble Constant parameter, Hq, represents the 
current value of the expansion rate for the universe. The 
k in Eq. ([l]) represents the three types of curvature for 
the spacetime of the universe as open, flat, or closed 
{k = —1, 0, 1, respectively). See Fig. [5] This is an 
inherent curvature of the empty spacetime itself, devoid 
of any matter or energy. However, as can be seen, the 
model does not use k directly, but rather the fractional 
curvature density parameter, flk- The total energy den- 
sity of the universe is split up into fractional pieces to 
represent its different compositional quantities for mat- 
ter, dark energy and curvature. The total energy den- 
sity, r^O: as measured today (z = 0), accounts for the 
sum total of all of the matter and energy in the universe. 
A critically dense universal model (typically accepted in 
cosmologjPH) , indicating that all matter and energy are 
accounted for, can be described by the relation 

no = n„, + nA + nk^i. (2) 

Note the curvature parameter. A:, has a negative sign 
(— ) originating from the GR spacetime equations of 
motiortlSi. However, the fractional curvature parameter, 
ilfe, has the relationship!^ 

oc -k, (3) 

such that 

{—1 : rife > Open (Negative Curvature) 
: = Flat (Zero Curvature) 

1 : fifc < Closed (Positive Curvature). 

The fractional matter density, £7™, represents the to- 
tal matter density in the universe. It can be separated 
into its constituents, flm = + ^c, into fractional 
baryon density, fib, and fractional cold dark matter den- 
sity, r^c, when the observation can constrain the distinc- 
tion (only BAG and CMB observation^. Note, for all 
values of ftm that were used in Sec. |llj we kept the con- 
stituents the same proportionate percentages of flm as 
the currently accepted values listed above. The A rep- 
resents the cosmological constant, the simplest model of 
dark energy. For changing the model of dark energy, 
the model uses an equation of state parameterized as 
w{z) = wq + Wa[z/{1 + z)pSl. The equation of state for 
dark energy includes, wq, which is the value measured 




FIG. 5: (color online) Universal Geometries (assuming 
isotropy and homogeneity) for illustrating various inherent 
spacetime geometries given by Qq (f^fc). TOP: Closed, i.e. 
sum of angles in triangle is greater than 180°; MIDDLE: 
Open, i.e. sum of angles in triangle is less than 180°; BOT- 
TOM: Flat, i.e. sum of angles in triangle is exactly 180°. 



today (z = 0), and its derivative, Wa, that allows dark 
energy to evolve in time, z, as the universe evolves. For 
the special condition when Wa = 0.0, the equation of 
state is constant, i.e. the density of dark energy does 
not change with time. Additionally, when wq = —1, this 
corresponds to the cosmological constant model. 



B. Experimental Observations 

Within the CosmoEJS package, we include 18 differ- 
ent experimental datasets for five different types of mea- 
surements: Type la supernovae (SNe), gamma-ray bursts 
(GRB), measurements of the Hubble expansion rate pa- 
rameter, H{z), baryon acoustic oscillations (BAG), and 
the cosmic microwave background (CMB) radiation. In 
cosmology, the evolution of the universe is modeled on 
scales too large to measure the evolution of a single 
galaxy or galaxy cluster from its formation to the present. 
Instead, cosmology combines observations of different 
galaxies — and the phenomena contained therein — at dif- 
ferent times and distances in their evolution to piece to- 
gether the dynamics of the universe. These different ob- 
servations and surveys provide independent and compli- 
mentary measurements of the expansion history of the 
universe and its composition. 

According to the big bang theory, the universe began 
from an initial state of extremely high temperature, den- 
sity and energy. When the universe had expanded and 
cooled enough for the photons to decouple from the pri- 
mordial soup of energy, they no longer scattered and were 
free to propagate throughout the universe, allowing for 
their detection. This surface is as far back as scientists 
can currently make measurements because all the mea- 
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surements involve some frequency of light. These high 
energy photons have been stretched with the expansion 
of the universe into microwaves. So, now the tempera- 
ture of the universe has cooled to about « 3 K. Working 
backwards, this gives a temperature of about « 380, 000 
K for the last scattering surface. We use the three fitting 
parameters for amplitude, and locations of the acous- 
tic peaks of the CMB temperature power spectrurrP: i) 
the acoustic scale, la (2*), 2) the shift parameter, 
3) the redshift of the surface of last scattering (SLS) of 
the CMB, Physically the acoustic scale and shift pa- 
rameter are the size, shape and position of the acoustic 
peaks of the CMB power spectrum for different values of 
the SLS redshift. The size, shape and positions of these 
acoustic peaks are very powerful in determining the val- 
ues of the cosmological parameters, due to the complexity 
of the many peaks in the CMB power spectrum and their 
ratios to each other. 

The SNe are standard "candles" (similar luminosity) 
used to form a redshift-distance relation, (distance 
modulus), to measure the rate of the expansion of the 
universe. According to these measurements, galaxies at 
large distances, in which the SNe reside, are receding less 
rapidly than Bubble's law would predict. This trans- 
lates to a slower expansion rate in the past, and that 
the nearby, later time SNe are expanding faster than the 
more distant, older SNe. Therefore, we are observing 
an accelerating universe. GRBs are added to fill the 
large void of redshift between the high-z SNe and the 
redshift of the CMB's last scattering surface, z « 1089. 
SNe are subject to extinction from the dust of the inter- 
stellar medium, however, GRBs are much brighter and, 
due to the high energy of gamma-ray photons, are rarely 
affected by the dust. The measurement of the GRB ex- 
tends the redshift-distance relationship to higher redshift, 
although, there is a redshift range of overlapping mea- 
surements for comparison and consistency. 

The measurements of the Hubble Parameter, ll(z) are 
an independent measurement of the expansion history of 
the universe. All of the other cosmological observations 
given in this article require an integration, but li{z) is an 
exact evaluation and comparison of Eq. ([T]) to the exper- 
imental data. In fact, Hi^z) is actually a direct measure- 
ment of the the differential age of the universe, Az/Ai, 
in other words measuring how the age of the universe 
changes as the redshift changes using the age differences 
of old elliptical galaxies that are passively evolvin^l^. 
They are used as standard "clocks" to directly probe the 
Hubble parameter. 

We compare the ratio of the sound horizon at the drag 
epoch, Tsizs), or when the baryons decoupled from the 
primordial universe to its effective distance, Dy(z) in 
the galaxy redshift surveys. This decoupling occurs at 
a slightly later time and lower redshift than the pho- 
ton decoupling because the baryons get "stuck" in grav- 
itational potential wells. The correlations in the galaxy 
redshift surveys consistently have a "bump" correspond- 
ing to the standard "ruler" measurement of the BAO. 



This measures the expansion of the primordial sound 
horizon in the galaxy redshift surveys. As an exam- 
ple of this physical ratio, the sound horizon at decou- 
pling in the range of ~ 153.19 Mpc, and effective 
distance D^{z = 0.57) = 2026.49 Mpc, for a ratio of 
rs/Dy — 0.076,'i^. Physically, this reflects the size and 
shape of the acoustic peak, and how it has evolved with 
the expansion of the universe. The BAO is specially 
suited for constraining {fif,, fid with galaxy clusters be- 
cause of the sensitivity of the sound horizon redshift to 
these parameters. 

Equation ([l]) can be used to compare several dark en- 
ergy cosmological models to these observations by allow- 
ing different values for the parameters, {f^m (= + 
ilc), ^^A, ^k,HQ, wq, and Wa}- While a model may flt 
one observation, cosmology involves the entire evolution 
of the universe, so it is important to use all the cosmo- 
logical observations, covering several redshift and cosmo- 
logical epochs, to find a best-fit model. 



IV. CONCLUSIONS AND FUTURE WORK 

CosmoEJS is a powerful new tool for cosmology edu- 
cation, and it is also precise enough to perform research 
grade calculations for testing most cosmological dark en- 
ergy models. They also allow the user to select inputs for 
parameters that are perhaps not scientifically accepted. 
This allows the user to discover how parameters influ- 
ence the shape of the curve for a particular theoretical 
model, thereby understanding the physical interpretation 
of a model's flt to the data. Variations of the programs 
have been used for science outreach and for classroom 
illustration. 

Future versions of the programs will involve a mini- 
mization method for the fitting of the cosmological mod- 
els to the data sets to provide best-fit cosmological pa- 
rameters. However, this will involve a different fitting 
method for each survey. We decided not to provide mini- 
mization in this version because this would distract from 
the pedagogical value of the program. Namely, when try- 
ing to find a best-fit model with a minimization routine, 
the user is not required to understand the physical inter- 
pretation of one fit over another. The fit is obtained by 
statistically comparing one model fit to another. Also, 
if the statistical fit is biased in some way, as explained 
in the examples above for fits, then the best fit could 
have unphysical parameter constraints. In the current 
version of the simulation, we are more concerned with 
understanding the physical interpretation of fitting par- 
ticular cosmological models to data sets. 

CosmoEJS allows non-specialists to manipulate cos- 
molgical models via their parameters and learn how to 
fit the model to experimental data sets. This manual pro- 
cess of changing parameter values also allows the user to 
see what parameter values do not fit the data. The pro- 
grams are useful for not only learning about cosmology 
but also data fitting itself, both visually and numerically. 
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The programs will continue to receive updates and mod- 
ifications for new, more precise data sets as these be- 
come publicly available. Using the ACDM model with 
this version of CosmoEJS, we find excellent fits to all 
the data sets with {Hq = 70.0 km/(s Mpc), = 0.045, 
= 0.225, = 0.27, Qa = 0.73, flk = 0.0, wq = -1.0, 

Wa = 0.0}. 



Appendix A: Facilitation with Easy Java 
Simulations (EJS) 

In order to test cosmological models, CosmoEJS is 
built with Easy Java Simulations (EJS)^. EJS is a Java- 
based software package for building interactive modeling 
simulations. The motivation for using EJS is that it com- 
bines a high performance coding language (Java) with 
easy graphical user interfaces (GUIs) and real-time plot- 
ting as the program calculates. The cosmological models 
for the evolution of the universe contain various parame- 
ters which relate the universe's physical composition and 
dynamics. CosmoEJS allows a user to select values for 
the model parameters and compare their personally con- 
structed model of the universe to the observational data. 
The programs allow the user to see the actual experi- 
mental data points and visually compare these data with 
the theoretical curves. The fitting of models with Cos- 
moEJS allows immediate feedback on a model's fit and 
physical interpretation both visually and quantitatively. 
Other methods of modeling in cosmology can require sev- 
eral programs and coding languages to produce resultJ^. 
Also, physical understanding can easly be missed when 
trying to learn complexities of several new programs. 
CosmoEJS is accurate and precise enough to test many 
cosmological models, however, its focus is to compare 
a few specific, physically relevant models that the user 
chooses. Unlike older programs, CosmoEJS allows the 
user to upload data sets, manipulate parameters, and 
plot several theoretical models; all in the same package. 
In addition, these programs allow the user to discover the 
cosmological model that best fits the current cosmologi- 
cal data sets. 

EJS is a populaili^, open source authoring tool which 
allows one to create Java simulations with minimal pro- 
gramming. A main advantage of EJS is that it allows the 
programs to be distributed along two fronts. 1) A simu- 
lation is packaged as a platform-independent executable 
".jar" file that can run on any system which has Java 
installed. There is no need for downloading compilers 
and learning new program languages in order to use this 
program. 2) The source code can be downloaded, and (as 
long as EJS is installed) the user can modify the code for 
their own particular usage, e.g. changing the integration 
method. CosmoEJS requires numerical integration us- 
ing the internal class. Integral, containing trapezoidal, 
Simpson's and Romberg's methods of approximating in- 
tegrals. We recognize that we are not using exact solu- 
tions, which may be different from these numerical ap- 



proximations. However , ex act solutions only exist for 
some very special case^^^^. We do provide the user the 
opportunity to change the method of integration and see 
convergence by adjusting the number of partitions and 
tolerance. Simply type "Integral . " and a list of meth- 
ods with proper input formats appears. See Fig. |6] 



O Description 6 Model O View 

OVariables C Initialization C Evolution C Fixed relations @ Custom O Elements 



{^■HMThe Friedmann model 


Integration 


Age 


BAD Primary 


CMBPrime 


Method to 


read data ► ' 





/// /Friedmann equations 

////Proper time Hubble in redshift H(z) 

public class InverseH implements Function { 

// This is uhere we define the model of the Universe 
public double evaIuate(double x) { 
// Unconunent the following lines when allowing 
// a varying dark energy 
double we = (-3 . 0*Ma*x) / ( 1 . 0+x) ; 
double wz - 3 . 0* ( 1 . 0+wO+wa) j 

// Also, uncomment the following additional 
// definition for qz when varying dark energy 
double qz = Math .pow( ( 1 . 0+x) , wz ) '^Math . exp (we] ; 

// Note: when wO=-l and wa=0, then qz = 1, 
// this is the cosmological constant model. 

//Friedmann model 

H = Hath.sqrt( (Omega_m*(1.0+x)* (1.0+x)* (1.0+x) ) 
+ Omega_v*qz + ((1 + x)*(l + x)* Omega_k)); 

return 1.0/H; 




M 

m 



O Description 6 Model O View 

O Variables C Initialization C Evolution C Fixed relations @ Custom O Elements 

' HWThe Friedmann model I nteg ration Age BAD Primarv CMBPrime Method to read d 



{///Routines for luminosity distance 
'public double integrate () { 

double dist; 

Function f - new InverseH (); 

// This is the integration of the model, now a function, 
//over the redshift to get the luminosity distance 

// Check each possible case of Curvature (Omegak) . 
if[Omega_k == 0) // If curvature is zero 

dist ^ hc^ (1 .0+z] ''Integral .romberg(f, zmin, z, partitions, tol); 
else^ // Other wise 

double OmegaKroot = Math . sqrt (Math . abs (Omegak) ) ; 

if [Omegak > 0) //If positive curvature 

dist - he* ( 1 .0^-z) * (1 .0/OmegaK_root) *Math . sinh{OmegaK_root« 
Integral . romberg ( f , zmin, z , partitions , tol ) ) ; 
else // If negative curvature 

dist = he* ( 1 .0+z) * (1 .O/OmegaKroot) *Hath . sin (OmegaKroot* 
Integral. romberg(f, zmin, z, partitions, tol)); 



} 

return 
} 



dist; //distance luminosity 



CDmmeni InlEgratior^s includE iLminositv distance, ADD, sound horizon 



FIG. 6: (color online) TOP: Screen shot of internal code 
showing the EJS methods used for numerically integrating the 
mverse Hubble expansion rate, Eq. BOTTOM: screen 

shot of internal code showing the EJS methods used for com- 
puting the luminosity distance for each situation of the piece- 
wise function using the Romberg's method. Note, in EJS it 
is not necessary for one to write a new code to perform the 
integration, nor is it necessary to include another file, since 
the Romberg routine is already included in the EJS numerics 
library. It is only necessary to call the integrand and sup- 
ply necessary input (integrand function, lower bound, upper 
bound, partitioning, and tolerance). 



Perhaps the most useful feature EJS offers (which dra- 
matically simplified this project) is the "View" window 
which provides users the means by which to easily and 
aesthetically construct GUIs, for input and output. A 
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user is able to "click and drag" a given feature and then 
manipulate its position to their preference, rather than 
taking the time to write the specific code for these fea- 
tures. For more on the general usage of EJS as both an 
executable " . j ar" file and the software package to build 
simulations, see Ref. ESI 



the other observational comparables on the successive ob- 
servational plotting tabs, the Hubble expansion rate, Eq. 
^ is used directly, H{z), and as the integrand in other 
observationfP. All integrations use the afore-mentioned 
Romberg's method already contained in EJS numerics li- 
brary, with appropriate partitions and tolerance, which 
the user can adjust to see convergence. 



Appendix B: Implementation of Cosmology 

In this appendix, we describe in more detail how the 
programs are designed to run and how cosmological data 
and theoretical equations have been implemented into the 
programs. For illustratative purposes, we begin with the 
implementation of the SNe data and theoretical equa- 
tions in the program containing all five observations^. 
The basic process by which a program runs consists of 
the program 1) loading observational data, then 2) us- 
ing parameter values to calculate and plot theorectical 
observables, and finally 3) comparing the fit with numer- 
ical value (values) to quantify the goodness of fit. 



a. Loading observational data 

Observational data of SNe are read in terms of red- 
shift value and distance modulus with experimental er- 
ror to set up the initial experimental data to fit to a 
model. Next, the program reads other observational 
data of GRB, BAO, H{z), and CM^. The most re- 
cently loaded data set is displayed on each observational 
data tab. There is also an option to supply a data set 
from the user^** . There are buttons to see the values of 
the experimental data loaded, and the theoretical values 
plotted on the screen after the program has calculated 
them^ 



b. Calculating the model and plotting 



c. Comparing the fit 

A visual comparison of a model's fit is evaluated by 
analyzing the plot and understanding why data points at 
certain redshifts match the theoretical curve and others 
do not. (See Fig. [l]-[3]for examples of models and their 
visual fit interpretations.) For numerical comparison, we 



Asne 



^grb 



^bao 



^crnb 



provide x , i.e. x — 

or likelihood, L = I\XiCo'w''^(xiXj)/\Xj with I\xi — 
x*-^ — xf'^ (if the inverse covariance matrix, Coy"^ {xiXj), 
is available), where C ~ We choose the -y? fit be- 

cause of its universal application with many different ob- 
servations including possible user supplied data set^. 
The is evaluated upon pressing the "CALCULATE! 
(Alt-C)" button, and is presented to the user before the 
plot finishes because this can be calculated needing only 
the theoretical comparable redshift values, rather than 
the actual plotting of all the theoretical data points re- 
quired for a smooth continuous plot. Messages display 
the x^ on the plotting panel on each tab for the partic- 
ular observation. A summary text output of the model 
parameter's values, CMB theoretical with actual exper- 
imental comparables, and total x^ or L is provided on 
the CMB tab (see Fig. [?]). 

Commit Microwave Background and Surrmarv 



Current model: omega_o = 1.000 critically Dense' (Fiat universe). 
HO = 70.000 km/s/Mpc H_0 = -1.000 M_a = 0.000 

Omegam = 0.270 Omegab = 0.045 Omegac = 0.225 

Omegak = 0.000 OmegaLambda = 0.730 

Age of Universe = 13.860 Gyr 



The user chooses initial values for the cosmological pa- 
rameter values and presses the "CALCULATE! (Alt-C)" 
button when ready to run the prograrrPSl. We plot a 
smooth continuous curve traversing all possible redshift, 
z, values. We use the evolution tab to step through all 
z values necessitated by the observations. The current 
z value is displayed in the "view" panel, with the same 
name. Given user inputs for the parameter values, the 
program generates a calculation of the Hubble expan- 
sion rate as a function of redshift, Eq. ([T]), (see Fig. |6]), 
and then performs numerical integration over all redshift 
values for the elliptical integral defining the luminosity 
distance, in the case of SNe and GRB, (see Fig. [6]). In 
order to compare to observations, the distance modulus 
is evaluated, and plotted versus redshift. To speed up the 
plotting, use the "spd" (steps per display) field, which re- 
flects how many steps between displayed frames and can 
be changed while the calculation is being performed. For 



Acoustic £cale [la) = 

Shift Parameter (R) = 

surface of Last scattering (z*) 



Theoretical 
304.672 
1.732 
1091 .42S 



Actual 
302.09 
1 .725 
1090.3 



SNe X"2 (X sne) = 562. S79 GRB X"2 (X grb) = 41.955 

H(z) X"2 (X_hz) = 21.901 BAO X"2 (X_bao) = 9.291 

CMB Likelihood (L) = 3.123 

Total Likelihood (L + X_sne/2 + X_grb/2 + X_hz/2 + X_bao/2) = 63S 



FIG. 7: (color online) The CMB output panel summarizes 
all the results, (color online matches the model), as well as, 
comparing the model to the CMB last scattering surface data. 
Here we show a model that fits all of the observations very 
well. 



Additionally, CosmoEJS allows users to deviate from 
critical density (fio = 1) which would be used to model 
a fiat universe. A message is displayed to notify the 
user whether the chosen parameters yield a open, fiat, 
or closed universe (see Figs. [iHs]). For critically dense 
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models, which are more commonly accepted, we incor- 
porate both a "derivation" and "fixation" option, with 
message. The derivation option allows the user to easily 
determine the value of one fractional density parameter 
based on the values of the other two. The fixation option 
holds a select parameter constant and allows the user to 
vary the ratio between the remaining parameters, such 
that their summation is always the critical density value 
minus the value assigned to the fixed parameter. Several 
models can be plotted on the same panel, each with its 
own color. Also, there are options to clear the plot of 
models, data, or everything, with the reset plot, data, or 
all button, respectively. For record keeping of the plots. 



the "Save Snapshot" button can be used for saving and 
printing^^ . 
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L. Barbato, Science 25 November 2011: pp. 1077-1078. 
See supplementary material at [URL will be inserted by 
AIP] for a table of all data sets included in this version of 
CosmoEJS with desriptions, as well as proper formatting 
for user supplied data sets. 

See supplementary material at [URL will be inserted by 
AIP] for a table describing the usage of all features in- 
cluded in this version of CosmoEJS. 
^® The redshift, z, is defined as 2: = {Xob3 — \emit) / Kmit, 
where \oba and Aemit are the observed and emitted wave- 
lengths, respectively, for the source''^. 
This is only comparatively speaking, for a precise calcu- 
lation of time from redshift, a model must be chosen, see 
Sec. IITBI 

The relationship between time and redshift (through a{t)) 
can only be calculated once a specific model is chosen. 
Actually, we evolve the scale factor, a, and calculate the 
time, t (in Gyrs), but switch the x and y axis, since this is 
the more common convention. 

It is important to compare the model both numerically and 
visually inspecting the plots, but also make sure the model 
is physical (e.g. Q.m = is unphysical, due to no matter). 
A detailed derivation of these models and an analysis for 
some special cases was recently made availableP^l. 
In cosmology, the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) spacetime metric!^ describes a homogenous and 
isotropic spacetime. 

When using the FLRW metric, the equation of motion is 
referred to as the Friedm ann equ ation. 
Recent data shows f^o ~ iHSHill^ although a more precise 
survejl^ uiay uncover a slight offset which can be modified 
in the source code, see Sec. Ill 

The use of the does have limitations. Namely, it can 
give biased results when used with unbalanced distribu- 
tions and errors of data. Specifically, lower redshift SNe 
dominate the fits due to a larger quantity of data points 
with smaller error bars. 



